High Fidelity Single Qubit Operations using Pulsed EPR 
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Systematic errors in spin rotation operations using simple RF pulses place severe limitations on 
the usefulness of the pulsed magnetic resonance methods in quantum computing applications. In 
particular, the fidelity of quantum logic operations performed on electron spin qubits falls well below 
the threshold for the application of quantum algorithms. Using three independent techniques, we 
demonstrate the use of composite pulses to improve this fidelity by several orders of magnitude. 
The observed high-fidelity operations are limited by pulse phase errors, but nevertheless fall within 
the limits required for the application of quantum error correction. 
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Pulsed magnetic resonance, which provides a way of 
manipulating quantum systems such as nuclear and elec- 
tron spins, has proved to be a powerful tool in the devel- 
opment of quantum computation. This is borne out by 
the early success of nuclear magnetic resonance (NMR) 
quantum computers, which have performed the largest- 
scale quantum computations to date [1, 2]. Always-on 
dipolar or exchange interactions are exploited to yield 
multi-qubit gates, whilst single-qubit operations are per- 
formed using classical radio frequency (RF) pulses. The 
scalability limitations surrounding NMR implementa- 
tions [3] that arise from the small nuclear Zeeman energy 
can be overcome by turning to electron paramagnetic res- 
onance (EPR) , analogous in many ways to NMR but with 
the advantage that pure states are experimentally acces- 
sible. 

These strengths have prompted many EPR-based solid 
state quantum information processing (QIP) propos- 
als [4-7]. The merit of such schemes is often argued on 
the basis of the decoherence time T2. However, the fi- 
delity with which operations can be performed also im- 
poses severe limitations on the viability of such proposals. 
We recently developed a methodology for characterising 
systematic errors in pulsed EPR, and used it to mea- 
sure typical errors in a commercial EPR spectrometer [8] . 
The most significant error present in a single-qubit op- 
eration is in the rotation angle, arising from spatial in- 
homogeneity in the pulsed RF field. This systematic er- 
ror in rotation angle is likely to persist even in the case 
of single-molecule EPR experiments, caused by miscali- 
brated control equipment. Fortunately, a number of ap- 
proaches to tackling different classes of systematic errors 
have been developed in the art of NMR, employing com- 
posite rotation sequences [9, 10]. A small subset of such 
approaches correct rotation operators rather than final 
states, and are therefore successful regardless of the ini- 
tial spin state. These approaches are applicable to quan- 
tum computing [11-14]. Of these, the BB1 sequence [11] 



exploits the precision in pulse phase control to correct for 
systematic errors in rotation angle, and is therefore ideal 
for our purposes. 

High-fidelity pulses are also beneficial to more tradi- 
tional EPR characterisation techniques such as correla- 
tion spectroscopy (e.g. 2D-HYSCORE [15]) by eliminat- 
ing or suppressing any spurious cross-peak signals and 
thus simplifying the spectral analysis. 

In this Letter we show how composite pulses can be 
applied in pulsed EPR to perform high-fidelity opera- 
tions on electron spins. We use these to demonstrate 
non-decaying Rabi oscillations, and provide further ev- 
idence using an error-sensitive electron spin echo enve- 
lope modulation (ESEEM) effect [16]. Finally, we use 
the error-measuring sequences described in [8] to esti- 
mate the residual error in the composite pulse. We find 
that in current EPR spectrometers, the chief limitation 
lies in the ability to accurately set the pulse phase, which 
restricts the effectiveness of the composite pulse. 

The paramagnetic species used here to perform er- 
ror measurements in pulsed EPR is i-NCgo (also known 
as N@C6o), consisting of an isolated nitrogen atom in 
the 4 S 3 / 2 electronic state incarcerated by a Ceo fullerene 
cage. It is an ideal system for these measurements be- 
cause of its extremely narrow EPR lincwidth and long 
relaxation time in liquid solution [17, 18]. T2 has been 
measured to be 80 /is at room temperature, rising to 
240 fis at 170 K [16]. 

The production and subsequent purification of i- 
NC60 is described elsewhere [19]. High-purity i-NCgo 
powder was dissolved in CS2 to a final concentration 
of 10 15 /cm 3 , freeze-pumped to remove oxygen, and fi- 
nally sealed in a quartz EPR tube. Samples were 0.7- 
1.4 cm long, and contained approximately 5 • 10 13 i- 
NC60 molecules. Pulsed EPR measurements were made 
at 190 K using an X-band Brukcr Elexsys580e spectrom- 
eter, equipped with a nitrogen-flow cryostat. 

i-NCeo has electron spin S = 3/2 coupled to the 14 N 



2 



nuclear spin 1=1. The EPR spectrum consists of three 
lines centered at electron g- factor g = 2.003 and split 
by a 14 N isotropic hyperfine interaction a = 0.56 mT in 
CS2 [20]. Most of the pulsed EPR experiments discussed 
below were performed using selective pulses on the cen- 
tral hyperfine line in the EPR triplet, corresponding to 
14 N nuclear spin projection Mj = 0, for which we can 
use a vector representation in visualising the evolution 
under RF pulses [8, 15]. The ESEEM experiment was 
performed on the Mj = — 1 hyperfine line, for which a 
full spin density matrix treatment is necessary [16]. 

In order to measure the quality of an operation inde- 
pendently of the starting state, we define the fidelity, J 7 , 
which compares the operator for the actual rotation with 
that of the ideal rotation. The fidelity, 



T= l -Tr(AB- 1 ) 



(1) 



takes a value between and 1 depending on how well the 
composite rotation B approximates the ideal rotation A 
(where A and B are unitary matrix operators). 

A general rotation of desired angle 9 with systematic 
error e, about an in-plane axis is given by 
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where a x and a y represent the Pauli spin operators. 
The BB1 corrective sequence has the form: 

Ho[e(l + e)] ft^Xl + e)] ft 02 [2^(l + e)] H^[n(l + e)]. 

(3) 

The fidelity of this composite pulse can then be expanded 
in e; all orders of the expansion up to and including co- 
efficients of e 5 equal zero for 
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For example, to achieve high-fidelity tt rotations, we 
choose 4>i = 0.580tt (104.5°) and 2 = 1.741tt (313.4°). 
Setting phases accurately is more difficult in EPR than 
in NMR, owing to the higher frequency; the fidelity of a 
BB1 7T pulse depends on small imperfections in the pulse 
phases, 50i and 502, approximately as 
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(O.75<50? - 1.12550i50 2 + O.5(50|) e 2 7r 2 , 
- (0.121501 - O.O9150 2 ) e 4 7T 4 + 0(e 6 ). 



(5) 



Qubit rotation is achieved in EPR through an on- 
resonance microwave pulse of controlled power and dura- 
tion. Rotation angle errors therefore arise from either 
pulse duration errors (which can be assumed uniform 
throughout the sample), or errors in the magnitude of 
the microwave field, B\, which varies across the sample 
depending on the homogeneity of the EPR cavity mode. 
For single- molecule manipulation there is no inhomogene- 
ity, but limitations in the resolution of power would still 
lead to systematic errors. 
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FIG. 1: Rabi oscillations for i-NC 6 o in CS 2 at 190 K (solid 
curve). BBl-Rabi oscillations exploiting BB1 composite 
pulses to remove the decay caused by pulsed field inhomo- 
geneity (dashed curve). 



The effect of B\ inhomogeneity is observed in the de- 
cay of Rabi oscillations as the RF pulse duration increases 
(see Fig. 1, solid curve). The inhomogeneity in B\ causes 
an inhomogeneity in the Rabi period, so spins in differ- 
ent parts of the sample gradually lose coherence under 
the influence of the RF pulse. By removing the error 
accumulated over these long rotations, BB1 composite 
pulses can be used to obtain Rabi oscillations that would 
be observed in the absence of an inhomogeneous RF field. 

For every desired rotation angle there are two dis- 
tinct phases which are required in the correction sequence 
(Eqn. 3). In order to be able to use the same two phases 
over the course of the experiment, the long RF pulses 
were divided into separate high-fidelity tt pulses, with 
the remainder provided by a simple pulse of length 9, 
where < 9 < tt. The BBl-Rabi oscillations shown in 
the dashed curve in Fig. 1 were obtained with the pulse 
sequence 

Tlo [0] ~ ( TloWi - "Ko. 587Tb 1 "] - ^1.747r[27r] - "K .58tt M )„ i 

(6) 

and demonstrate that this sequence does not accumu- 
late errors arising from a B\ inhomogeneity. Since T2 ~ 
200 /is, which is long compared with one Rabi time pe- 
riod, electron spin decoherence is small on the timescale 
shown in Fig. 1. 

In a recently reported mechanism for ESEEM, the fre- 
quency content of the echo modulation was found to have 
a clear dependence on the fidelity of the echo-refocusing 
pulse [16] . The echo decay is modulated by two frequency 
components, 5 and 25, whose relative amplitudes are a 
function of the rotation angle of the refocusing pulse. For 
a 7r(l + e) = tt + 9 e refocusing pulse, where 9 e represents 
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FIG. 2: Fourier transform of ESEEM in i-NC6o shows two 
frequency components 26 and 52 kHz. For simple tt refocusing 
pulses, the 26 kHz component is attributed to pulse error 
(solid line) and is removed when a BB1 composite tt pulse is 
used (dashed line). For a refocusing pulse of approximately 
0.617T, the 52 kHz component, which would be present due to 
pulse error, can also be removed by using a BB1 composite 
0.617T pulse (dotted line). 



the (small) absolute error in the refocusing pulse, the 
relative magnitudes of the two components is given by 



F(2S) 



= 20?. 
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However, for a refocusing pulse of twice the magic an- 
gle [9]: 



6 = 2 cos" 

the S component dominates: 

F(S) _V2 
F(2S) 6 e ' 



(8) 
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Fig. 2 compares the Fourier Transform of the echo decay 
using either simple or BB1 composite pulses, for nominal 
refocusing pulses 9 — tt and 6 = 0.6087T. For the latter ro- 
tation angle, phases <f>\ = 0.5497T and <p2 = 1.6467T were 
used in the BB1 pulse. When the error-compensated 
pulse is applied, the secondary frequency component is 
removed in each case, further demonstrating the ability 
of this composite pulse to correct for systematic errors 
in rotation angle. It is worth noting that in these ES- 
EEM experiments, the BB1 "refocusing" pulses operate 
on spins which are dispersed in the rotating x-y plane. 
This illustrates the effectiveness of the BB1 composite 
pulse over a range of initial states. 

Rotation angle errors can be measured by compar- 
ing the rates of echo decay in two multi-pulse se- 
quences: Carr-Purcell (CP) and Carr-Purcell-Meiboom- 
Gill (CPMG) [8, 9]. The echo decay in CPMG shows no 
sensitivity to rotation angle errors (after every even num- 
bered cycle) and is dictated purely by decoherence. Echo 
magnitudes in CP show a cumulative sensitivity to errors 
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FIG. 3: Comparison of the echo signal decays in the reference 
CPMG sequence (dots) with that of the error-sensitive CP 
(crosses) sequence provides a measure of rotation angle errors. 
CP echoes generated with simple n pulses decay more rapidly 
owing to rotation angle errors. When the simple tt pulses 
are replaced with BB1 composite n pulses (open circles), the 
decay rate is very close to that in CPMG. 



in the refocusing tt pulses and hence decay faster than in 
CPMG, as shown in Fig. 3. By performing a comparison 
as described in Ref . [8] , we estimated that the refocusing 
pulses had a standard error of 0.17T (18°), equivalent to 
a fidelity of T = 0.988. 

By replacing the tt pulses in CP with BB1 composite 
tt pulses we are able to make an estimate of the residual 
rotation angle error in the high-fidelity rotation. The 
error accumulated in the CP sequence is dependent on 
the number of cycles applied and hence the sensitivity of 
this technique is limited only by the maximum number 
of pulses that can be applied in one experiment (with the 
spectrometer used for these experiments, this is currently 
32, which implies fewer than 8 BB1 composite rotations). 
From the decay of the BB1 corrected CP sequence, we 
are able to conclude that the rotation angle error is at 
most 0.017T (2°), corresponding to a fidelity T > 0.9993). 
The phases used in this BB1 experiment were measured 
using a SPAM sequence [8] and found to be <\>\ = (0.587± 
0.008)tt and (f> 2 = (1.742 ± 0.011)tt. Using Eqn. 5 these 
imply an expected fidelity T = 0.9999. 

Error correcting codes will be crucial to the operation 
of a quantum computer [21], compensating for both de- 
coherence errors and the gate operation errors described 
here. The threshold of error probability for these codes 
varies as a function of the overhead, but generally lies be- 
tween 10" 3 and 10" 4 , or a fidelity greater than 0.999 [22]. 
Despite the inherent 10% systematic error, the compos- 
ite pulses described here are capable of producing opera- 
tions which meet this threshold fidelity. However, errors 
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in pulse phase cause the fidelity to fall short of the the- 
oretical optimum of 1 — 1CP 6 . In addition to the BB1 
composite pulse described here, longer composite pulses 
have been proposed which can arbitrarily reduce the sen- 
sitivity to systematic error [13]. In order for these higher- 
order composite pulses to be effective in pulsed EPR, 
greater phase control is required. An analysis of random 
errors present would complement the current work in de- 
termining where the compromise between lengthy pulse 
sequences and error reduction lies. 
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